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PREFACE 


The contents of this document are designed to serve two purposes. First, 
to provide automated computer Identification of minerals and compounds from 
unknown samples for those staff members (mineralogists and geochemists) In the 
Planetary Science and Applications Branch, Ames Research Center, who have been 
using the more traditional manual Identification technique employing the Fink 
Index. Second, to provide detailed Instructions and worked examples for use 
In graduate level courses In mineralogy (Geol 211) and X-ray analysis applica- 
tions (Geol 255) at San Jose State University. This document presumes the 
reader has either had prior course work or direct experience in mineral and 
rock type Identification employing X-ray diffraction and emission techniques. 

Appendix I, a compendium of computer printouts describing worked examples 
discussed In Section IV of this report. Is available from either Dr. Marshall 
Haddock (San Jose State University, Department of Geology) or Maxwell Blanchard 
(Ames Research Center, Planetary Science and Applications Branch). 
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IDENTIFYING UNKNOWN MINERALS AND COMPOUNDS FROM X-RAY 
DIFFMCTION PATTERNS USING THE JOHNSON AND 
VAND FORTRAN IV COMPUTER PROGRAM 
Frank T. Kyte* and Maxwell B. Blanchard 
Ames Research Center 


I. INTRODUCTION 


The Johnson and Vand FORTRAN 1\ Program for Identification of Multiphase 
Powder Diffraction Data can be a useful tool to interpret complex diffraction 
patterns from mixtures of minerals. The program compares d-spacings and inten- 
sities of unknown diffraction patterns with every compound in the ASTM X-ray 
diffraction file. It then prints out every compound that satisfies the input 
parameters and data. These parameters Include; chemical information, inten- 
sity match requirements, allowable error in pattern measurement, number and 26 
range of d-spacings observed, and background intensity. The 10 standards 
having the best match are then subtracted from the original unknown pattern, 
leaving a residual pattern. 

This report describes the program operation and how to interpret the out- 
put. Input parameters are discussed in enough detail that a new user can 
understand why patterns are accepted or rejected by the computer. A principal 
reference for this report is "FORTRAN IV Program (Version 10) for the Identifi- 
cation of Multiphase Unknown Powder Diffraction Patterns" by G. G. Johnson, Jr. 
and V. Vand, published and distributed by ASTM, 1969. 


II. JOHNSON - VAND PROGRAM DESCRIPTION 


A. Input Pp’-ameters 

Successful use of this program depends on the choice of input parameters. 
The most important of these are the error window and the log intensity match, 
which define how closely a standard pattern must match the unknown in its 
d-spacings and intensities. The parameters log intensity match, number, and 
percent of line matches apply only to the range of d-spacings within the input 
high and low range of 20. These five parameters interact in a very compli- 
cated way. When the d-spacing range is unusually narrow, number and percent 
should probably be decreased and log intensity match increased. 

(1) Subfile option- The subfiles are alloys, organics, inorganics, miner- 
als, and dual (inorganic compounds). The subfile option allows the operator 
to limit the search, and thus the output, to subfiles in which he is interested. 


*Geology Department, San Jose State University, San Jose, Calif. 95192. 


(2) Resolution- This parameter instructs the program that a high- 
resolution Guriier camera has been used and the intensities are "fuzzed" to give 
a "normal" resolution. This parameter is not used for Debeye-Scherrer diffrac- 
tion patterns or for diffractometer scans. 

(2) Error windoi}- The error window is defined by a variable d* which, 
unlike d(X), has a linear relationship with 0 and best approximates a con- 
stant error window for the observed range of 9 values. The error window is 
a very important parameter because it defines how closf ly a standard line must 
match a d-spacing in the unknown pattern. 

Errors in the diffracted line measurement (in Angstroms) are roughly pro- 
portional to d^. A constant error window is best approximated by transforma- 
tion to reciprocal spacings that are 1/d. When used by the program, the 
reciprocal spacings are multiplied by 1000, rounded down to the nearest inte- 
ger, and added to 0.5. Thus, input d-spacings from 200 X to 0.6 A are repre- 
sented by a set of numbers from 5.5 to 1666.5 (e.g., 5.5, 6.5, 7.5, . . ., 
1666.5). This new number, d*, allows the program to perform Integer arithmetic 
for much faster calculations and defines a constant e’’ror window for any value 
of 0. 

d* = l/d(A) X 1000 

Figure 1 shows a plot of d(X) and 1* versus for CuKa = 1.5418. 


ai 



d* 


Figure 1.- Plot of d* and dA versus °0 for CuKa = 1.5418. (To convert 
from dA to d* read d in A from upper scale using line labeled dA 
which corresponds to a °0 value on left hand scale. Using this °0 value 
read the corresponding d* using the line labeled d* whose value is given 
on the lower scale. For example, a d-value of 2.98A corresponds to an angle 
of 15° 0 and is equivalent to 335 d*.) 
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To determine the optimum error window for the Instruments at Ames 
Research Center, a mixture of approximately equal weights of flourite, galena, 
rutile, chromite, and olivine was X-rayed using the dlf fractomet:»r and the 
Debeye-Scherrer camera. 

To see how much error could be expected, d* was calculated for multiple 
measurements on the Debeye-Scherrer powder films and for each of the diffrac- 
tometer scans shown in Tables 1 and 2, respectively. The diffractometer scans 
were corrected because the goniometer was consistently about 0.05“ 0 too high. 
For 99 percent (508 measurements) of the line measurements, the maximum varia- 
tion was 4d*, or an error of ±2d*. Of the remaining 1 percent (5 measurements) 
of line measurements, the maximum variation was 6d*. 

A nominal error window of ±2 would seem to be appropriate for most 
identifications. However, inaccuracies in the ASTM file itself, or variations 
in d-spacings caused by solid solution, may necessitate the use of a larger 
window (■' } or gicater). One should be cautious of using too large an error 
window because it can produce a large number of possibilities, and, in some 
cases, the windows could even overlap, creating a continuous coverage of 
several degrees 20. . Jually, a broad error window of ±4 is suitable only for 

data in conju ction with chemical information. 


TABLE 1.- DEBEYE-SCHERRER CAMERA FILMS 


Sample 

pattern 

Camera 

radius 

Sample 
grain size 

Quality 
of lines 

Times 

measured 

Total 

lines 

Date 

ASTM 1 

114.6 mm 

coarse 

spotty 

4 

120 

7/22/75 

ASTM 2 

114.6 mm 

fine 

continuous 

4 

112 

7/24/75 

ASTM 2A 

57.3 mm 

fine 

continuous 

4 

66 

7/28/75 

ASTM 3 

57.3 mm 

coarse 

spotty 

4 

95 

7/29/75 

ASTM 5 

57.3 mm 

fine 

continuous 

1 

— 

8/19/75 


TABLE 2.- DIFFRACTOMETER SCANS 


Sample 
pat tern 

“20/min 

cps 

' 

Slits 

Peaks 

Date 

ASTM 

4A 

1 

250 

1/2 

14 

7/30/75 

ASTM 

B 

1 

100 

1/2 

18 

8/4/75 

ASTM 

C 

1 

250 

1 

22 

8/4/75 

ASTM 

D 

1/2 

250 

1 

22 

8/4/75 

ASTM 

F 

1/2 

250 

1/2 

13 

8/5/75 

ASTM 

F 

1/4 

100 

1/2 

31 

8/5/75 

BSTM 

1 

1/4 

250 

1 

— 

10/3/75 

BSTM 

2 

1 

250 

.J_.j 

— 

10/9/75 
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(4) Log intenaity match- This very important parameter Is a requirement 
for Intensity agreement between the unknown and an ASTM standard. All inten- 
sities are input on a 3-dlgit scale of 1 to 100, but all program calculations 
are on a 1-dlglt log scale of 0 to 9. The 3-digit (I3) Intensities are con- 
verted to the 1-diglt di) log scale by the formula 


ll - 5 log 1 3 


The conversion scale thus becomes 


I3 

0-1 

2 3 4-6 

7-10 

11-15 

16-25 

26-39 

40-63 64-100 

II 

0 

12 3 

4 

5 

6 

7 

8 9 


Thus, an input intensity of AO is equivalent to an input of 63 and very differ- 
ent from an input of 39. Program users should keep this in mind when inputting 
intensities. 


Log intensity natch is defined as 

log intensity match (1 1 ) “ 9 - (I ^STD ~ ^2 UNK.)mAX 

where (I^STD “ ^2 UNK)mAX the maximum intensity disagreement between any 
line in the standard and its counterpart in the unknown pattern. When select- 
ing an input value for log intensity match, the user must remember that any 
standard pattern with a match of this value or less will be rejected as a 
possibility. 

For example, when log intensity match is input with a value of 2(13) it 
is converted to a log scale value of IClj). and then intensities of every 
line in the standard are compared with corresponding intensities of lines in 
the unknown (which may be zero if no line is present in the unknown to compare 
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with one in the standard pattern) , and log intensity match is calculated for 
each line in the standard. If the log intensity match of any single line in 
the pattern has a value of l(Ii) or lower, then the standard will be rejected 
as a possibility. Therefore, in this example, it is necessary that 


^ ■ ^^ISTD " ^lUNK^MAX ^ ^ 

or that 

^^ISTD " ^lUNK^MAX 8 or 9 

Essentially what this means, in returning to the I 3 scale, is that lines 
with intensities of 40 to 100 in the standard must match with corresponding 
lines in the .inknown. Failure to do so results in rejection, even if every 
other line ir the standard pattern except one, say with an intensity of 40, 
matches perlectly. 

(5) Nwr.h.tr- This is the minimum number of line matches that the program 
requires for in identification. If this minimum is set too high, a pattern 
may not be if'entified (e.g., diamond has only 5 lines). 

(6) Percent of line matches- This parameter allows the operator to select 
the minimum percent of line matches between an unknown and a standard. In 
most cases, 33 percent should be sufficient, requiring that at least 1/3 of 
the lines in the standard must match the unknown pattern. 

(7) Icycle- This parameter instructs the program to go through a second 
cycle, during which time, all previous parameters are Increased. The amount 

of increase is specified by another set of parameters described in Section III. 
In general, it is suggested that this parameter not be used because all param- 
eters must be increased. If Icycle were used to Increase the error window, a 
larger number of standard lines could be compared, but rejection of previous 
matches could occur from stricter requirements for log intensity match or 
number. 

(8) Observed ranae of d-spacings- This is defined by two parameters, DHI 

and DLO, which are the range of d-spacings capable of being observed in the 
diffraction experiment. DHI must be higher than the highest input d-spacing. 
For example, a 3A line might be the largest observed spacing when the experi- 
mental range was up to 20X. If DHI were Input as rather than 20X, then 

standards with strong lines between 3 and 20A could be matched when they 
should have been rejected. DLO should equal the lowest d-spaclng in the 
unknown pattern. Note that back reflections (from Debeye-Scherrer camera 
patterns) are often omitted in the AbTM file. It is desirable to Include all 
of the unknown patterns down to about I.IA. 

(9) Background intensity- This parameter instructs the program to dis- 
regard all lines in the standard pattern of this value or smaller. This is 
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difficult to quantify on powder films and should be given as about half the 
lowest Intensity observed. On diffractometer scans, it is a good idea to set 
the counter scale so that the most intense peak(s) is at about 100 on the 
chart. Then the background (as well as other peaks) can be read directly off 
the chart. 

(20) Option- This parameter uses chemistry. If the input is MAJOR, all 
the compounds (regardless of chemistry) will be printed. If the input is CHEM 
or SOLID, then only compounds which are chemically correct will be printed in 
the output as potential identifications. If the input is MINOR: (a) all com- 

pounds fitting parameters for a match will be printed, and (b) if a compound 
is chemically correct, the intensity agreement will not be applied. 

(22) Chemical information- Chemical information can be used in this 
program. It is possible to list all elements and functional groups which are 
known to be present and those which are known t*" be absent. Chemical informa- 
tion is highly desirable because it may narrow down a list of over 100 possi- 
bilities to less than 10 or 20. It is strongly suggested that chemistry be 
used whenever possible. 

(22) Preconceived possibilities- This option is to be used for phases 
which are definitely known to be present. It can be very useful in that part 
of the output (described below) where standard phases are subtracted from the 
unknown pattern. Preconceived possibilities are always subtracted first. 

Thus, if a match was selected on the basis of certain lines in the unknown 
pattern generated by the preconceived possibility, when it is subtracted from 
the residual unknown pattern, it will have fewer lines to compare with and it 
will be scaled down to have a minimal effect on the residual pattern. This 
will be discussed in more detail in the next section. 


B. Program Output 

The program output comprises of five sections: INPUT, OUTPUT, REPORT, 

MATCHD, and SUBSTI. Figures 2 through 8 show a portion of the output obtained 
by Johnson and Vand as their own example. It was generated by a known mixture 
of equal weights of Cap 2 , Si02, ZnS, CaC 03 , and PbS. The pattern was obtained 
from a Norelco diffractometer at l®/min, using Cu radiation and peak 
intensities . 

(2) INPUT- This is a summary and description of the input parameters and 
data (fig. 2), and should always be checked carefully for correctness. If 
this section is not checked, the user might miss an error comment such as, 

"YOU MAY HAVE KEYPUNCHED INPUT DATA WRONG - A SPACING IS OUT OF ORDER - 
PROGRAM ORDERS AND WILL CONTINUE." 

(2) OUTPUT- The OUTPUT section (fig. 3) lists every compound, regardless 
of chemical information (unless CHEM or SOLID options are used), that fits the 
input data and parameters. It also lists the various numerical measures of 
degree of fit oi the unknown to the standard patterns, as well as a reliabil- 
ity factor which is a composite index calculated from the percentage of line 
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for any possible error messages.) 
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Figure 3.- OUTPUT. (This section lists all the possibilities which fit the input 
parameters. All Intensities are I 3 values.) 
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Figure 3.- Concluded. 
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Figure A.- Continued. 



l.b4« 

6- 

3 1.JJ9 

t- 

3 

l.tlt 

C— 

3 u.ass 

6- 

3 

i.C2* 

O- 

J C.S39 

t- 


C.S77 

i- 

^ O.aZi 


« 

I.C^e 

2— 

< C»cJl 

3- 


l.l»i 

i- 

j l.C«9 

w 

d 

1. 18C 


i C.90S 

i- 

1 

l.«< 1 

i- 

< C.iSl 

i- 

C 

1.^09 

3- 

2 



<.t9! 

3- 

2 



d.9«* 

^- 

1 



i.iU 

2- 

1 



I.Ci« 


1 





1 





1 



^ I.J40 


1 



i.9l9 


i 



1. 


1 




C.Sfc 1- c 

C.Se» 1- 0 

C.S.C I- C 

1- C 

l.Coi 1- O 

1. 1^4 1- b 

1.14< 1- C 

i- 0 

l.iki 1 - c 

1.2i< I- C 


Figure 4 



6- 

3 1.775 

6- 

1.432 

6“ 

3 1.903 

6- 

1 

6- 

3 2.245 

6- 

1.787 

6- 

i 3.063 

6- 

1.871 

6“ 

3 1.951 

6 - 

2.181 

6- 

i 5.831 

6- 

2.442 

6- 

3 2. 195 

2- 

2 .998 

6- 

3 


6.844 

6- 

3 


1.9v9 

2- 

1 


1.6 79 

2- 

1 


1.708 

2- 

1 


1.896 

2- 

1 


1.932 

2- 

1 


2.C9S 

2- 

1 


2. 681 

2- 

1 



Concluded 


c 


ixe l-tLlOhlKC TAdCF MtPHtS l;»T^ P«R<LLEL UP the UMUNOMN PATTEKN UITh EhE BESI FUUNl. IN IHE MEPOHT 

NU ATItNPI l> NAOE lU SLU IMAt. I IMEKSlIlt! 

iNE «STP pupbep thL Tpt peliaeilitv eactcn eead The Column up each stanoaho 


Ire PKECUnCtlVtU PCAildlLI IlcS HAVE A AEllAtlCIIV FACICH OF 944 


IFE 

SIAAIIPC 

kITF 

AklP 

PUPCEf 

5C59< 

IS 

Ft 

S 


IFE 

5IANUANU 

ktfn 

ASIP 

ncpbei* 

5C45E 

I* 

SI 

C2 


iFE 

5IANIARC 

■ I7f 

AkIP 

NUH BeH 

4(E«4 

IS 

CA 

F2 


|MC 

SIAKJArU 

kllM 

AkIP 

kUPBEH 

3C560 

IS 

IS 

S 


TFE 

5IANUAMU 

H ITH 

ASIN 

NCMMEk 

SC!EC 

1 S 

CA 

C CJ 


iFE 

SIAKCAftC 

■ nr 

ASTP 

kUMBEM 

IdCAiA 

IS 

IH 

S 


IHE 

ilAKUARU 

bllH 

ASIP 

kCPBEP 

5C4S2 

IS 

IP 

s 


IFE 

5 IANCARC 

M IIF 

ASIR 

ncmdem 

I2CI4I 

1 s 

S 



IF.E 

5IAACAMC 

kITF 

ASIP 

PUPBEA 

5056 5 

IS 

SI 


IFE 

S lANOAMd 

kllH 

ASIP 

PCPdeH 

16(1(1 

1 s 

FE2 

!S1 S4 



kicei) IlCdSi CIEFFCAi I ICMICNEMS 

UIHII OllCI 

AEIN ASIP 

^0!9^ »049(l 

75 4J 


A .E45 
6.EA5 
S.831 

A.41i 

4.^/i.C 14 A.idll 4^<di 

4.;ie4 J9-^5 

i.<«C 

J.BOCC b A.OHi5 J.Bibl 

J.ES4 

J.75S 

).7itt 

J.431 lUd-ui 


ASIM 

ASIP 

ASIM 

ASTM 

ASTM 

ASTM 

ASTM 

ASTM 

kOdbA 

50566 

5C586 

100434 

50492 

130141 

50565 

IBOTOO 

42 

35 

35 

32 

29 

16 

13 

13 


i5-l5 
4 - 2 

o- i 

6- 2 


15- Id 

15-10 

15-10 

25-15 


Figure 5.— MATCHD. (The 10 best matches are listed In order of decreasing reliability 
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pattern.) 
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Figure 7.** Concluded. 
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Figure 8.- SUBSTI FINAL REPORT. (The 10 best possibilities ere listed with their seeling 
factors. All the subtractions of SUBSTI are listed in parallel.) 
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Figure 8.** Concluded 


matches. Intensity agreement, and spacing agreement. Experience has shovm 
that a reliability factor greater than 33 and a log Intensity match greater 
than 25 Indicate a strong probability of the presence of that phase. If the 
input parameters are too restrictive (e.g., error window of 2 and log Intensity 
match of 4 when the pattern has only 6 to 8 lines) and no possibilities match 
the Input data, then an error message will be printed in this section. It can 
be Identified by the statement "ABEND CODE IS: SYSTEM 00C5." 

(Z) REPORT- The entire pattern for those compounds which are chemically 
correct is given (fig. 4). If no chemical Information Is provided, then the 
10 compounds with the highest reliability factors will be printed. Each com- 
pound listed is headed by its ASTM number, reliability factor, formula, and to 

the right is a number denoting Its subfile, as follows: 

blank » An inorganic pattern 

2 * Organic 

3 “ Dual 

4 = Inorganic, mineral 

5 “ Greater than loA (organic) 

6 = Organic, mineral 

7 = Greater than lOX (Inorganic) 

8 = Cross-Indexed minerals 

9 = Greater than lOX (dual) 

If chemical information is used as input but no chemical formula occurs 
on the ASTM card, then the compound will not be listed in REPORT, MATCHD, or 
SUBSTI. Thus, the OUTPUT section should be visually scanned for possibilities 
with high reliability factors and intensity matches, but without chemical 
formulae. 

(4) MATCHD- The 10 best matches are listed in older of decreasing relia- 
bility (fig. 5) . The numbers and asterisk to the rig'uc denote the subfile and 
data quality. Below this to the left, the unknown pattern, with intensities 
and error windows, is listed. The next column lists the d-spacings of the 10 
best matches which are then listed to the right. Each possibility is headed 
by its ASTM number and reliability factor. Under the heading are the 1 3 
intensities (grouped according to the Ij scale) which correspond to the 
standard d-spacings to their left. 

(5) SUBSTI- This section attempts to subtract standard phases from the 
unknown pattern, one at a time, based on intensities. This results in a final 
residual pattern. Intensities of the unknown are converted to the same log 
scale as the standards (figs. 6 and 7) before execution. SUBSTI FINAL REPORT 
(fig. 8 ) lists d-spacings and scaled intensities of all the standard together. 

A scaling factor is used to scale down the intensities of the standards 
to those of the unknown. This is calculated as 

.... ^ . 2 ^unknown 

Scaling factor ■ ■= 

2 ^ ^standard 
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where lunknovn intensities of the unknown pattern whose d-spacings 

match those of the standard to which it is compared. These values may be nega- 
tive if already subtracted by a previous standard with a higher degree of con- 
fidence. It is here that a preconceived possibility is useful. A known phase 
will be subtracted first, scaling down any later possible phases to Insignifi- 
cance if their identification relied on lines belonging to the known phase. 
SUBSTI FINAL REPORT has the same format as MATCHD, except that now scaling fac- 
tors have been calculated, standard intensities have been scaled down, and each 
standard has been subtracted from the unknown pattern. 


III. PROGRAM INPUTS AND OPERATION 


The Johnson-Vand program is simple to operate. A minimum of 8 data cards 
(described below) are placed inside a deck of 15 job cards. The data follow 
the first 14 job cards and the final job card is placed at the back of the 
deck. The deck is then given a job header card (see fig. 9). Table 3 
describes the data card formats. 

TABLE 3.- DATA CARD FORMATS 


DATA CARD 1 
DATA CARD 2 


Col. 1(A) 


Col. 3(1) 


This card may contain any information (not beyond column 76) 
to identify the run. 

The parameter card. The program will supply a reasonable 
set of parameters, default parameters if a value of zero or 
blank is given. The parameters are; 

IFL, the subfile of the ASTM powder diffraction file to be 
searched. Set IFL as follows: 

blank * all subfiles 
A » alloy 

0 = organic 

1 = inorganic 

D = dual (e.g., both Inorganic and organic) 
default = blank (entire file) 

M = mineral 

RESOLt resolution condition 

0 » normal 

1 = Guinier 
default = 0 


Col. 5-8(1) MM, the number of lines in unknown pattern 

default = computer will count number of peaks 


24 


TABLE 3.- Continued. 


Col. 9-12(1) 

IWy the error window 
default ■ 2 

Col. 13-16(1) 

KUT0FFt the minimum acceptable value of leg intensity match 
default ■ 2 

Col. 17-20(1) 

NUMBER^ the minimum acceptable number of line matches 
default “ 3 

Col. 25-28(1) 

IPERCT, the minimum acceptable number of line matches 
default ■ 33 

Col. 29-32(1) 

ICYCLE, the number of cycles 
default * 0 

Col. 33-36(1) 

IWCY, the increase of IW per cycle 
default ■ +1 

Col. 37-40(1) 

KVTOCY t the increase of KUT0FF per cycle 
default = +5 

Col. 41-44(1) 

NUMCYf the Increase of NUMBER per cycle 
default = +1 

Col. 45-48(1) 

IPERCY, the increase of IPERCT per cycle 
default ■ +10 

Col. 49-52(1) 

JBKGR, the background intensity (0 to 100 scale) 
default = +1 or half lowest given intensity 

Col. 53-60(F) 

DHI, the upper limit of the d-spacing range observed in the 
diffraction pattern 

default “ error window above largest measured d 

Col. 61-68(F) 

DLOt the lower limit of the d-spacing range observed 
default * error window below smallest measured d 
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TABLE 3.- Continued. 


Col. 69 


DATA CARD 


DATA CARD 


DATA CARD 


73(A) OPTION 

SOLID = Solid solution 
MINOR = Minor phase 
MAJOR “ Major component 

default * Major 

3 The positive element and preconceived possibilities card. 

Chemical elementa on this card and on Card 4 must be right 
justified, two columns per element; that is, single letter 
elements must be punched with a blank, then the letter: 
for example, a list of sodium, hydrogen, lithium, iron, 
nickel, oxygen, and platinum would be punched 

NA HLIFENI OPT 

With the first N in column 1 and none beyond column 50. 

RE (Rare Earth) punched in columns 49 and 50 will cause all 
the rare earth elements from La to Lu, inclusive, to be 
used as positive elements, in addition to any others punched 
in columns 1 through 48. The letters TU (transuranium) 
punched in columns 49 and 50 of the position element card 
will cause all the elements from Ac to Fm to be included 
automatically in the positive element list. A Q in column 1 
calls for the periodic table. 

Preconceived possibilities- The ASTM card numbers patterns 
which are to be printed out even though not found by the 
search program. Each must consist of six digits with no 
spaces. (A maximum of five patterns can be listed as 
preconceived possibilities. As an example, if patterns 1-24, 
6-432, 11-4, 12-100, and 16-432, were known to be present, 
then the card would contain 010024060432110004120100160432 
(starting in column 51). In other words, the six numbers 
are divided into a two-digit field for the card set and a 
four-digit field for the card in that set. The minimum 
contents for this card is 000000.) 

4 The negative element card. These must be punched right 
justified, two columns per element as in Card 3 with none 
beyond column 80. A Q punched in column 1 will cause all 
elements (through Fermium) not specifically designated as 
positive elements to be used as negative elements. 

5 Positive Functional Groups. Those functional groups which 
are known to be present in the unknown (from spot or other 
chemical tests). 
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TABLE 3.“ Concluded. 


DATA CARD 6 


DATA CARD 7 
thru N 


Negative Functional Groups. Those functional groups which 
are known to be absent in the unknown (from spot or other 
chemical tests) . 

NOTE: The input for both these cards are as follows: 

Punch chemical group exactly as the group appears in 
figure 3 (***OUTPUT***) under FORMULA. Notice the spaces - 
this is extremely Important. Eight groups can be external, 
starting in columns 1, 11, 21, etc. 

These cards contain the pairs of d-spacings and intensities 
measured on the unknown pattern. These need not be in 
order. Punch the d-spacing, then its intensity (scale 
1-100). The d-spacing mist contain a decimal point and the 
Intensity must not contain a decimal point. These cards are 
of the format (8(F7 . 3, 13) ) . Then one blank card. Cards 1-N 
may be repeated as many times as unknown diffraction 
patterns are available. 


A printout can be obtained with only a data set of spacings and intensi- 
ties rcard 7-N) along with 7 blank cards in the proper places. The program 
provides a reasonable set of parameters (by default) if none are input. 

On Card 2, all numerical parameters should be right justified in their 
fields. For example, an error window (IW) of 3 which has a field of col- 
umns 9-12 should be punched into column 12. Remember that a zero on this card 
is read as a blank and a default parameter will be used. It is impossible to 
give a parameter a value of zero. 

On Cards 7-N, the diffraction data cards, the format is (8(F7.3,I3)) . 

That is, for each card which has 80 columns, there are 8 fields of 10 columns 
each for paired spacings and intensities. Each of these 10 column fields is 
divided into a 7- and a 3-column field for the spacing and intensity, respec- 
tively. In the 7-column field, if no decimal is used, all blanks are counted 
as zeros and a decimal is placed automatically before the three right-hand 
digits. If a decimal is used, the spacing value may be placed anywhere within 
this field. The 3-column intensity field must be occupied by an integer which 
must be right justified. All blanks in this field are read as zeros. At 
least one of the intensities must equal 100. 

All data cards 1-N must be followed by a blank card to end the search. 

A typical run would probably use the mineral subfile. The error window 
would be defaulted to 2, or set at 3. All other parameters would probably be 
defaulted except perhaps the log Intensity match, which might be set at 4 (or 
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Figure 9.- Illustration of input cards. 
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Figure 9.- Continued. 
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7 for a good diffractometer scan). If an identification did not result, data 
should be rechecked and the error window should be exp. ided and other [ aram- 
eters increased by small amounts to limit the amount of output. With experi- 
ence, an operator should get a feel for what the capabilities of the system 
are. 


IV. A WORKED EXAMPLE 


A set of worked examples has been prepared. The reader is referred to 
Appendix I which is a package of computer printouts. This package is avail- 
able on loan from Max Blanchard at Ames Research Center or Dr. Marshall E. 
Haddock at San Jose State University. Appendix I comprises nine printouts 
obtained from Johnson and Vand's example (figs. 2-8) and from two mixtures of 
minerals (Tables 1 and 2) which were X-rayed with the Debeyc-Scherrer powder 
cameras (57.3- and 114.7-mm diameters) and the X-ray diffractometer (at 
different scan speeds). 


A. Sample Number 5 (CaF 2 , SIO 2 , XnS, CaC 03 , PbS 

This is the entire output of the Johnson and Vand example. The entire 
Inorganic file was searched. Chemical information was useu. The only input 
parameters were the range of observed d-spacings. The effect of using chemi- 
cal information can easily bo appreciated here as it reduced a total of 170 
possibii . ties to only II. From there, it easily Identified the five phases 
presto.; . 


B. ASTM 2, 2A, 5, and 4F 

These four printouts were obtained from a mixture of approximately equal 
amounts of PbS, Cap 2 , Ti 02 , (Mg, Fe) ^SiOi, , and Fe(Cr,Al) 204 . Unlike the 
Johnson and Vand example, two minerals In this mixture are subject to consid- 
erable solid solution and standard patterns in the Powder Diffraction File, for 
these minerals represent only points along a solid-solution sequence. Thus, 
standard intensities and d-spaclngs could vary significantly from those 
observed in an unknown sample. 

Data for ASTM 2 were obtained from a large Debeye-Scherrer powder camera. 
All of the parameters wore defaulted and the mineral subfile was searched. 

Cap 2 and PbS were easily identified In OUTPUT with high reliabilities and 
intensity matches. TiOo had the next highest reliability but a low-intensity 
nuitch. (Mg,Fe) 2 Si 04 and Fe (Cr , A1 ) 2 O 4 were not identified because their strong 
lines at 2.2?A and 2.0b8X, respectively, were not measured. Inspection of the 
X-ray film confirmed their existence but their intensities were very weak and 
spotty. The data were not good enough to identify these phases. In SUBSTl 
FINAL REPORT, cordierite and montmorl 1 Ionite were listed as good possibilities, 
as well as the other three ph.ases. Here, chemical information would have been 
helpful in determining the actual phases present. The program was, however. 
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successful, using default parameters, In searching the approximately 2000 
mineral subfile and listed five good possibilities without the use of chemical 
Information. Three of these five were correct. 

Data for ASTM 2A were obtained from a short exposure on a small Debeye- 
Scherrer camera. All the parameters were defaulted and the mineral subfile 
was searched. None of the possibilities listed had confidence levels high 
enough to be considered as an Identification. It Is assumed that the exposure 
time was too short to darken the film and that too few diffraction lines were 
measured for an Identification. Thus another, longer, exposure was taken and 
named ASTM 5. 

Output for ASTM 5 was from data for 28 lines, whereas ASTM 2A had only 17 
lines for the same range of d-spaclngs. The data were run at an error window 
of ±3, and only the mineral subfile was searched. Results were comparable to 
those on the large camera. Ti02, PbS, and CaF 2 were easily Identified, but 
the two ccmponents with solid solution were missed again. Attempts at broaden- 
ing the error window were futile because too many possibilities were selected. 
In general, a broad window of ±4 is suitable only for data with chemical 
Information. 

Data for printout ASTM 4F was obtained from a slow diffractometer scan at 
1/4** 20/m1n. This was much more successful than with the powder cameras. This 
is attr: ted to the diffractometer's ability to resolve weak lines and the 

greater .. curacy of line location and peak height which are inherently better 
on diffraction charts. These data were processed with default parameters and 
only the mineral subfile was searched. Without chemical information only 20 
possibilities were listed in OUTPUT. Six of the 20 represented the 5 known 
mineral phases and 2 others (spinels) were isostructural with the chromite. 
SUBSTI FINAL REPORT lists all except olivine as having a high probability, and 
the only incorrect phase listed there is BeO, an unlikely component. 

Observation of the log- intensity matches listed led to the conclusion that 
the output could be improved if this parameter were increased. Another com- 
puter run was made with a log- intensity match of 4. This is ASTM 4FA. A 
brief glance at OUTPUT and SUBSTI FINAL REPORT indicates an excellent match of 
all phases except olivine, with the only other phases listed being BeO and 
NiFe 204 , again, a spinel isostructural with chromite. 


C. BSTM 1 and 2 

Another mixture, this time composed of equal weights of rutile (Ti02) , 
ga era (PbS), zircon (ZrSi 04 ), Flourite (CaFo), and calcite (CaC 03 ) was pre- 
parea. Of these, only the calcite exhibits solid solution with the replace- 
ment of Ca by Mg, and it was a relatively Ca-pure sample. Two diffractometer 
scans were run at 1/4° 20/min and at 1° 20/min and labeled BSTM 1 and BSTM 2, 
respectively. 

Both data sets were recorded in two ways. One was to read directly off 
the diffraction charts (BSTM 1 and 2) and the other (BSTM lA and 2A) was to 


subtract 0.05° 0 to compensate for misalignment of the goniometer. (Note: 

0.05° 0 Is equivalent to about Id*.) These four data sets were then Input ^ 
with default parameters except for log-intensity match at A and DHI at 8.8A5A 
corresponding to an angle of 10° 20. Background radiation was read directly 
off the charts and only the mineral subfile was searched. 

BSTM 1 at 1/4° 20/mln had an output of only 10 possibilities uncorrected 
and 13 when corrected by 0.05° 0. In both cases, all five known phases were 
Identified. The corrected output had higher reliability factors In every case 
except zircon, which matched with 11 lines Instead of 12 after being corrected. 

In both cases, SUBSTI FINAL REPORT listed the five known phases as the 
only ones with scaling factors above 0.1, the lowest of these being calcite at 
0.206. The remaining possibilities had been scaled down so low because their 
identifications had been dependent on lines generated by those phases with 
higher reliabilities. When it came time to calculate a scaling f'.ctor for 
these incorrect phases, there was little or no pattern left for them to 
compare to. 

BSTM 2 at 1° 20/mln uncorrected had an output of eight possibilities, 
correctly identifying every phase except calcite. The corrected data set 
listed only 10 possibilities, including calcite, and reliability factors for 
the known phases were greatly enhanced. Once again, in both cases, the known 
phases were the only ones listed with scaling factors greater than 0.1 and a 
low factor for calcite at 0.189. Some very successful runs indeed! 


IV. SUMMARY AND SUGGESTIONS 


It has been demonstrated that this program can easily identify mixtures 
of minerals which can be quite difficult to resolve manually. There is still 
some difficulty in identifying phases which are subject to solid solution, it 
is hoped that when the master file is updated from Set 19 to Set 24 this will 
be less of a problem. The update will nearly double the data file and this 
should be accomplished i" the near future. When a number of spacings have not 
been identified, solid solution should be suspected. The best way to resolve 
this would be to use ch» mical information, a wider error window (perhaps ±4), 
and the OPTION parameter. Remember, though, that it is possible that a phase 
is not identified because its pattern is not in the master file. 

When inputting data, default parameters should be sufficient in most 
cases. An error window of ±3 might be needed wh'^n there is a poor resolution 
of lines. With confidence in intensity measurements a log-intensity match of 
4 will greatly reduce the number of possibilities, requiring that standard 
intensities down to 26 must match with lines of the unknown pattern. A second 
cycle is usually not desirable! 

Remember to input the entire range of d-spacings, especially DHI. When 
inputting intensities, keep in mind the log scale the program uses. Note that, 
in the worked examples, all intensities are input on this scale. Remember, at 
least one intensity rmst equal 100. 
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As an operator uses the program more often, he will become more familiar 
with its capabilities. He should always keep in mind, however, that the pro- 
gram merely lists possibilities with different degrees of confidence. The 
worked examples presented were quite successful; however, at other times incor- 
rect identifications with high reliabilities have been observed. In the end, 
it is the responsibility of the operator, not the computer, to decide which 
phases are actually present. 
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